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In a preceding paper I lj , Mukhopadhyay and I studied the diffusive motion of a tagged molecule 
in an heterogeneous glass-forming liquid at temperatures just above a glass transition. Among other 
features of this system, we postulated a relation between heterogeneity and stretched-exponential 
decay of correlations, and we also confirmed that systems of this kind generally exhibit non-Gaussian 
diffusion on intermediate length and time scales. Here I extend this analysis to higher temperatures 
approaching the point where the heterogeneities disappear and thermal activation barriers become 
small. I start by modifying the continuous-time random- walk theory proposed in [lj], and supplement 
this analysis with an extension of the excitation-chain theory of glass dynamics. I also use a key result 
from the shear-transformation-zone theory of viscous deformation of amorphous materials. Elements 
of each of these theories are then used to interpret experimental data for ortho-terphenyl, specifially, 
the diffusion and viscosity coefficients and neutron scattering measurements of the self intermediate 
scattering function. Reconciling the theory with these data sets provides insights about the crossover 
between super-Arrhenius and Arrhenius dynamics, length scales of spatial heterogeneities, violation 
of the Stokes- Einstein relation in glass-forming liquids, and the origin of stretched exponential decay 
of correlations. 



I. INTRODUCTION 

Many of the most important questions in the physics of 
glassy materials have to do with spatial heterogeneities. 
What is the scale of these heterogeneities? In what sense 
are they dynamic as opposed to structural phenomena? 
What is their relation to stretched exponential decay of 
correlations? To anomalous, non-Gaussian diffusion? To 
violations of the Stokes-Einstein relation between diffu- 
sivity and viscosity? I address each of these issues in this 
paper by studying a relatively simple model in which all 
of them arise, and which allows me to make contact with 
experimental data. 

In a preceding paper [lj, Mukhopadhyay and I stud- 
ied a model of the diffusive motion of a tagged molecule 
in a heterogeneous glass-forming liquid. This model is 
based in part on the excitation-chain theory of the glass 
transition. [2, 3] It consists of two different kinds of molec- 
ular environments that we called "glassy" and "mobile." 
To a first approximation, the tagged molecule is either 
frozen at a fixed position in a glassy domain or is mov- 
ing diffusively in a mobile region. The glassy domains 
are compact regions of relatively low energy and entropy, 
surrounded by more highly disordered regions in which 
molecular rearrangements can occur. Because of the mo- 
bility of the molecules within and at the boundaries of the 
mobile regions, these boundaries diffuse slowly through- 
out the system on an "a" relaxation time scale, and the 
system as a whole is ergodic. A tagged molecule in a 
glassy domain is stationary until it is encountered by a 
moving boundary, at which time its motion is enabled 
and it diffuses until it finds itself back in a glassy do- 
main. In this sense, the model has elements of kinetically 
constrained systems, d i, H, 0, 0, i, El 

The analysis in [l[ pertained only to behavior at tem- 
peratures T low enough that the heterogeneities are much 



larger than the molecular spacing. With that assump- 
tion, we found a clear decoupling between the slow a 
relaxation and the faster mobile motions that we identi- 
fied as "/?" cage-breaking rearrangements. Fluctuations 
in the sizes of the glassy domains produced stretched- 
exponential decay of the self intermediate scattering 
function at long times, and the distribution over spatial 
displacements exhibited an anomalously broad exponen- 
tial (non-Gaussian) tail during the transition from 0- to 
a-like motion. 

My purpose here is to reformulate [l[ in a way that ex- 
plores the decrease in the size of the heterogeneities as a 
function of increasing temperature above the glass tran- 
sition. My ultimate goal is to understand the crossover 
between the thermally activated phenomena that occur 
near the glass temperature and the non- activated, fluid- 
like transport that is described at high temperatures by 
mode-coupling theory. 

I first reformulate the preceding analysis in SeclTIl 
pointing out how the theory must be modified as the sys- 
tem changes from being predominantly glasslike at low 
temperatures to predominantly mobile near the transi- 
tion from super-Arrhenius to Arrhenius dynamics at a 
crossover temperature Ta- Then, in Sec JIIH I review 
the excitation-chain theory 0, S] of these temperature- 
dependent quantities and show how, in a simple approx- 
imation, the basic ideas of that theory may be extended 
to describe the behavior near the crossover. 

The resulting temperature-dependent scattering func- 
tions and diffusion and viscosity coefficients are com- 
pared with experimental data for ortho-terphenyl (OTP) 
in Sec lIVI Among other results, this analysis makes 
it appear that the spatial heterogeneities in OTP are 
relatively small, and that the observed violation of the 
Stokes-Einstein relation near the glass temperature oc- 
curs because the mechanisms responsible for diffusion 
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and viscosity become different from one another as the 
glass-forming liquid becomes more like a solid. Section 
Ivl contains a discussion of the way in which stretched 
exponential decays of correlations emerge in a theory of 
this kind. Finally, in Scc lVIl I summarize my present 
opinions about the issues raised in the first paragraph of 
this Introduction. 



II. CONTINUOUS-TIME RANDOM WALKS 

The continuous-time random-walk (CTRW) analysis 
fill IT2I used in [l[ required that molecular motions con- 
sist of just two alternating modes - glassy and mobile - 
each with its own characteristic time scale. The present 
analysis is based on a similar two-stage assumption, mod- 
ified in a way that more naturally crosses over to high- 
temperature behavior in which the tagged molecule un- 
dergoes only normal diffusion. The two-stage assumption 
is a serious one because it precludes exploring molecu- 
lar motions in the continuous range of short time scales 
between molecular vibration periods and (3 relaxation, 
i.e. in the so-called "ballistic" regime. This analysis 
cannot resolve motions on those very short time scales, 
or on length scales appreciably smaller than the inter- 
molecular spacing. Nevertheless, the CTRW technique 
is uniquely capable of describing the larger scale motions 
which do exhibit a clear decoupling between fast and slow 
modes. The power of this kind of analysis is illustrated 
in [l3j |. which has served as a starting point for much of 
the present work. 

With this short-time limitation in mind, denote the 
two temperature-dependent time scales by T a (T) for the 
long, a relaxation time, and tm{T) for the ordinarily 
shorter, mobile relaxation time. As in [l|, let the char- 
acteristic size of the heterogeneities be £ R* (T) , where £ 
is a length of the order of the average molecular spac- 
ing and R* is a dimensionless function of T. Then the 
characteristic fluctuation time for the heterogenities is 

T«(T)=R* 2 (T)T a (T). (2.1) 

A central ingredient in this analysis is the probability 
Vg{T) that the tagged moleecule is in a glassy domain. 
This probability is a geometrical quantity, equal to the 
fraction of the system occupied by glassy domains, which 
is assumed to be approximately a constant in time al- 
though the positions of the domains themselves are fluc- 
tuating. 

The waiting-time probability distribution for a 
molecule initially in a glassy domain was found in [l[ 
to be 

^ G {t*)^2e- 2 ^, (2.2) 

where t* = t/r^ and t is the time in dimensional 
units. Note that ipcit*) 1S normalized over t* rather 
than t. This distribution, which has the Kohlrausch form 
exp (—const, x t* b ) with a stretched-exponential index 



b = 1/2, was obtained by assuming a Gaussian distribu- 
tion of domain sizes R with a scale size R* (T) . The index 
b is probably a temperature dependent quantity, rising to 
unity as T increases to the crossover temperature Ta] but 
here, as in [l[, I make no attempt to estimate this T de- 
pendence and simply keep b = 1/2, independent of T. 
As will be shown in Sec|V] this "intrinsic" value of b is 
generally not what is observed experimentally. 

The main departures from [l[ that are needed here 
concern molecular motions in the mobile regions. The 
revised model must make a smooth transition between 
the low and high temperature regimes. At low T, the 
analysis in 0, y] suggests that the system consists pri- 
marily of glassy domains surrounded by thin, tenuous 
mobile regions. When a molecule enters a mobile region 
under these conditions, I assume that it moves a distance 
proportional to R* and, with high probability, is back in 
a glassy domain at the end of this motion. By letting this 
jump-length scale with R* , the size of the heterogeneities, 
I account roughly for fact that the mobile regions may ex- 
tend over large distances, even though they occupy only 
a small fraction of the total volume at low temperatures. 

In contrast, at higher temperatures near Ta, the sys- 
tem consists primarily of mobile regions in which small 
glassy domains are embedded. In this case, I assume 
that a molecule makes multiple diffusive jumps in the 
mobile region, and only with low probability finds itself 
in a glassy domain after any single jump. A constraint 
on the model is that, as T approaches and exceeds Ta, 
molecular motion should cross over from anomalous to 
normal diffusion with a diffusion constant Dm oc £ 2 /tm- 

Because no stretched-exponential behavior is expected 
to be associated with diffusion in mobile regions, let the 
waiting time between mobile jumps have a simple expo- 
nential distribution: 

V>M(t*) = Ae- A '*, (2.3) 

where A is the dimensionless ratio between the glassy 
and mobile residence times. There are various, equally 
plausible choices for A, the simplest of which is 

A = t q /t m . (2.4) 

With this choice, and with the definition of t* , the charac- 
teristic mobile waiting time is of the order of tm R* 2 oc 
£ 2 R* 2 /Dm- If the time between jumps scales like the 
square of the domain size, (£R*) 2 , then the characteris- 
tic jump length must scale like £R*; therefore write the 
jump-size probability in the form 

W = WaW 2 6XP ' (2 - 5) 

where r* = r/(£ R*) is the jump length r in units of £ R*, 
a is a dimensionless parameter, and d is the dimension- 
ality. (Ordinarily d = 3 in the examples to be shown 
here. For some purposes such as in Sec. IIII1 however, it 
is useful to keep general d.) 
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This picture so far is consistent with the idea that, 
when R* is large at low T, the jump size and the wait- 
ing time are both large. At high T however, where R* 
is small, the picture makes sense only if multiple jumps 
can occur in mobile regions. Thus, assume that a mobile 
jump is followed by another mobile jump with probabil- 
ity Vm = 1 — V G) an d that a molecule ends in a glassy 
domain after such a jump with probability V G . This as- 
sumption means that, at low T where Vg is nearly unity 
and Vm is small, the molecule makes mostly infrequent, 
long jumps; whereas, at high T where Vm is close to unity 
and Vg is small, the molecule makes repeated, fast, short 
jumps. 

The mathematics of the CTRW analysis is much the 
same here as in [l[ . As usual, it is easiest to compute with 
the Fourier-Laplace transforms of the constituent prob- 
ability distributions. Define two different distribution 
functions, n G (r* ,t*) and %()■*,<*), for molecules start- 
ing, respectively, in glassy domains or mobile regions, and 
moving distances r* in times t* . Their Fourier-Laplace 
transforms are 

ftj(*,u) = J e- lk ' r 'dr* e- ut ' nj {r*,t*)dt\ (2.6) 

where j — G,M, k* = k^i?*, and k is the dimensional 
Fourier wave vector. 

The Laplace transform of the glassy waiting time dis- 
tribution ipc(t*) is 

i>a(u) = / e- ut 'i> G (t*)dt*. (2.7) 
Jo 

As discussed in this function has an essential singu- 
larity at u — and a branch cut conveniently chosen to 
lie along the negative real u axis. The probability that 
the molecule has not yet left its glassy domain at time t* 
is 

^ G (t*')dt*' = (1 + 2 VF)e- 2V ^; (2.8) 
and the Laplace transform of this function is 



cj> G {u) 



1 - i>G(u) 



(2.9) 



The analogous functions for mobile motion are 

^m(u) = A t ; 4>m(u) = 



u + A' 



u + A 



(2.10) 



The Fourier transform of the jump-length probability dis- 
tribution in Eq. (|2.5[) is 



/(**) = exp - 



2 L*2 



a 2 k 



(2.11) 



where k* — |k*|. 

The next step in the analysis is to consider a sequence 
of jumps in which a molecule starts in a mobile region 



and does not enter a glassy domain until its last jump. 
The corresponding propagator is 



T{k*,u) 



V G f(k*)$M(u) 
l-V M f{k*)^ M {u) 



(2.12) 



A related function describes the sequence in which the 
molecule remains in the mobile region after its last jump: 



Q{k*,u) 



4>m(u) 



l-V M f(k*)ip M (u) 



(2.13) 



With these ingredients, the functions nc(fc*,u) and 
fiM(k*,u) become 



h G (k*,u) 



and 



n M (k*,u) = 



4>g{u) + G(k* ,u)\jj G (u) 

1-T{k*,u) # a (u) 
1 N G (k*,u) _ 
u W(k*,u) ' 

4> G (u)F{k*,u) + g(k*,u) 

l-f(k*,u)i> G (u) 
1 N M (k*,u) _ 
u W(k*,u) ' 



(2.14) 



(2.15) 



where 



N G (k*,u)= \ i-v M f{k*)\ [i +«/ A ; 

(2.16) 

N M {k*,u) = [l-^ G (u)l V G f(k*)+u/A; (2.17) 



and 



W(k*, u) = 1 - V M f(k*) + u/A - V G f(k*) ^ G {u). 

(2.18) 

In Eqs. (|2.14j) and (12.15|) . the first term in the numerator 
corresponds to the case in which the molecule ends in a 
glassy domain and the second term to the case where it 
ends in a mobile region. 

The Laplace transformed self intermediate scattering 
function, F s (k*,u), is the weighted average of n G and 
n M - 

F s (k*,u)=V G (T)h G (k*,u)+V M {T)n M {k*,u). (2.19) 

The scattering function that is measured experimentally 
is the inverse Laplace transform: 



F s (k*,t*) 



-\-ioo 



du 
2?r7 



F g (k*,u). 



(2.20) 



As in this integration is performed by closing the 
contour around the branch cut on the negative real u 
axis. The crucial ingredient is, for u — > — w± ie, e — > +0, 

POO 

4> G (-w ± ie) = -4 / f d£ e - w ^ ±2ii 
Jo 
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w 



2y¥ 

7 ,.3/2 



-l/t, 



Erfi 



A(w) TiB(w), 



(2.21) 



where Erfi is the imaginary error function. The imagi- 
nary part of ^g{~w ± ie) is its discontinuity across the 
branch cut on the negative u axis. With this formula, it is 
straightforward to compute the discontinuity across the 
cut for the functions hj(k*,u), j = G, M in Eqs. (|2.14p 
and (|2.15[) . and in this way to compute the scattering 
functions by closing the contour around this cut. That 



hj(k*,t*) = lim 



dw ■ 



x Im 



t* 

7T W 

-w + ie) 



_ W{k*, -w + ie) 



j = G, M. (2.22) 



The integrands are well behaved asw^ +0, therefore it 
is possible to let the lower limit of integration be w = 0. 



III. EXCITATION-CHAIN THEORY 

Evaluation of the formulas in Sec|TT] requires explicit 
expressions for the temperature-dependent quantities 
r a (T) , t m (T), R*{T), and V G {T). These expressions 
are needed not just as asymptotic forms near the glass 
transition but, more importantly, as physically motivated 
approximations valid near the crossover from anomalous 
to normal diffusion near Ta- That crossover coincides 
with a transition from super- Arrhcnius to Arrhcnius dy- 
namics, and from heterogeneous to homogeneous spatial 
properties. 

I propose to use the excitation-chain (XC) theory 0, Q 
to derive these formulas. So far as I know, it is at present 
the only theory of glass-forming liquids that relates their 
dynamic and thermodynamic properties directly to the 
behavior of real molecules with short-range interactions, 
and that does so in a manner consistent with the ba- 
sic principles of statistical mechanics. This theory re- 
mains speculative and incomplete. It does not even tell 
us whether the glass transition is thermodynamically well 
defined or whether, as I strongly suspect, relaxation times 
grow increasingly rapidly but smoothly all the way down 
to zero temperature. If the latter is true, then the XC 
theory describes just the onset of that process and, in 
doing so, gives us just a clue about what energy scales 
and collective motions might be involved. In either case, 
it provides a framework for understanding what might 
be happening near Ta- 

In order to go beyond the immediate vicinity of the 
glass temperature, I need to extend the published XC 
results, and therefore need to provide a brief summary 
of the basic ingredients of this theory. The fundamen- 
tal question to be addressed is why molecular rearrange- 
ments that create or annihilate density fluctuations are 
not simply thermally activated Arrhenius processes. My 



answer is that simple, one-step, activated processes do 
occur, but that they do not always lead to stable, con- 
figurational rearrangements. Think of the creation of a 
density fluctuation very roughly as the formation of the 
glassy analog of a "vacancy-interstitial pair" in which a 
molecule moves out of an energetically favored site to 
a neighboring, less favored position, perhaps leaving a 
lower density at the original site and producing a higher 
density at the new site. Denote the energetic cost of 
this move by k B Tz ■ Then think about what happens 
during the next thermal fluctuation. The energetically 
favored next move is recombination; the interstitial falls 
back into the vacancy. At high enough temperatures, 
however, this energetic advantage is outweighed by the 
entropic advantage of there being many different next- 
neighbor sites available. The vacancy and the interstitial 
in effect dissociate from each other in a single activated 
step, and the resulting density fluctuation is an Arrhenius 
process with a temperature-independent energy barrier. 

At lower temperatures, by definition below Ta, disso- 
ciation must become a multi-step process. The simplest 
such mechanism that I can imagine is a chain-like se- 
quence in which one molecule pushes or pulls its neighbor, 
that neighbor displaces the next neighbor, etc., finally 
leaving a vacancy at one end of the chain of displacements 
and an interstitial at the other. Chainlike (or stringlikc) 
excitations appear to be ubiquitous in the neighborhood 
of jamming transitions in amorphous materials. For ex- 
ample, see the work of Glotzer et al. [HI, HH, [HI . The 
energy cost of such an excitation must, to a first approx- 
imation, be proportional to the length of the chain; but 
this cost may be compensated by an entropic advantage; 
the longer the chain of displacements, the more routes 
are available for it to take. 

To interpret this picture mathematically, compute the 
rate of thermal activation of an excitation chain consist- 
ing of N molecular steps and extending a distance R 
in units t, where R is dimensionless and i is the same 
molecular length scale that was introduced in Eq. (|2.ip . 
This problem is closely related to - but not exactly the 
same as - the problem of computing the probability of 
iV-step random walks with extensions R. Write this rate 
in the form of an attempt frequency, say t^ 1 , times an 
activation factor exp(— AG(N, R)/k B T), where: 



AG(N, R) = k B T z + Ae - k B T In IT + E VI 



(3.1) 



The first term on the right-hand side, k B Tz, is the bare 
activation energy discussed above. The remaining terms 
describe the statistics of the chain. However, I omit the 
diffusion term, proportional to R 2 /2N, from the usual 
random-walk analysis because it turns out to be negligi- 
ble for present purposes. 

In the second term in Eq. (|3.ip . eo, is the activation 
energy per link in the chain, roughly the energy required 
to move a pair of molecules far enough away from each 
other to allow a third molecule to pass between them. 

The third term, W(N, R), is a sum over chain config- 
urations weighted by the disorder of the glassy environ- 
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ment in which these excitations occur. It has the form: 

(3.2) 



In W(N, R) « vN-—^-R, 



where exp (y) is the number of choices that the successive 
links in the chain can make at each step. The function 
j{T) is the mean-square amplitude of the fluctuations in 

e /k B T; 



(3.3) 



To = eo/vkB is the characteristic temperature deter- 
mined by the energy eg. With the approximations used 
here, To turns out to be equal to both the Kauzmann and 
the Vogel-Fulcher temperatures in the XC theory. j /£ 
is the inverse localization length associated with diffu- 
sion in disordered systems. Because glassy disorder is on 
molecular scales, 70 seems likely to be about unity; but a 
first-principles derivation remains one of the most uncer- 
tain and speculative parts of this theory. The derivation 
given in almost certainly breaks down at large N and 
low temperatures. Thus, it seems likely that this calcu- 
lation makes sense only at temperatures not too close to 
To, and that even the existence of a nonzero Tq remains 
problematic. 

The last term in Eq. (|3.f p . Ei nt , is especially impor- 
tant for the present analysis. This is the energy that 
makes it unfavorable for the links of the chain to lie near 
each other. In @, Q, I used Flory's approximation [l7| 
for the self-exclusion energy of a polymer chain in a d- 
dimensional environment: 



(3.4) 



where the temperature T int sets the scale of this energy 
up to a dimensionless factor of the order of unity. This 
approximation is valid, however, only in the limit of large 
N and R; whereas we need a theory that makes sense in 
both the latter limit and for vanishing values of those 
quantities near T = T4. For the moment, therefore, 
rewrite this energy in the form 



E int (N, R) « k B T i7 



J(N) 
R d 



(3.5) 



where the as-yet undetermined function J(N) is equal to 
TV 2 in the limit of large N but becomes a much weaker 
repulsion in the limit N — *■ 0. 

To compute the activation rate from Eq. (|3.1[) . first find 
the value of R = R* at which AG (AT, R) is a minimum. 
The resulting function AG* (TV) = AG(N,R*) is conve- 
niently written as 



AG* (TV) 



T Z +T&(T, TV), 



(3.6) 



where 



a{T,N) = -vN 



T 



and 



+(<f+l) 



R* 



2d 



2 d Ti n t 

TTj(T)T 



Tin 



J(N) 



J(N) 



d+1 



(3.7) 



(3.8) 



If T > To, and if the second term on the right-hand 
side of Eq. (l3.7|) is sublinear in N, then a(T,N) has a 
maximum at TV = N*(T) which, in analogy to nucle- 
ation theory, is the length of the critically large chain 
that nucleates a long-lasting density fluctuation. When 
a chain fluctuates to a size larger than TV*, it is highly 
likely to continue growing without bound, and thus to 
dissociate the vacancy- interstitial pair. Thus the activa- 
tion energy for this process is AG(N* , R*), and the a 
relaxation time r a is given by 



In 



TQ 



Tz 
T 



a(T), 



(3.9) 



where a(T) = a[T,N*(T)}. 

For temperatures low enough that TV* ^> 1, the Flory 
limit in Eq. (|3.4p is accurate; that is, J(N) w TV 2 , so that 



R*(T) 



N*(T) 



nj(T)TT int ~ 


1 


2 


2 

d-1 


2d 









, (3.10) 



7T 7 (T)T 

2d 



d-1 1 
d-1 



v (T - To) 



d + 1 
d-1 



(3.11) 



and 



a(T) a (d-1) 



"tt7(T)" 


d 
d-1 


~2{TT lnt )f 


2d 




v{T-T ) 



2d 



R*(T). 



(3.12) 



Thus we recover the Vogel-Fulcher formula for r a (T) near 
T = T , with a(T) cx (T — Tq)" 1 for d = 3. 

Note that, in Eq. (|3.12p . a(T) is linearly proportional 
to R*(T) independent of dimensionality d. This formula 
and parts of the preceding analysis resemble recent re- 
sults of Eckmann and Procaccia [1 81 ] , who make a math- 
ematically systematic study of a two-dimensional model 
and find both chainlike excitations and a linear relation 
between the super- Arrhcnius activation energy and the 
spatial extent of the chains. 

Now return to Eq. (|3.7p and notice that, if J(N) ~ 
cN d+1 for small N, this equation can be made to have 
the form 



n(T. V ) = .V (1-^ 



(3.13) 
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by choosing 



(2d) d [v(T A -T )] 



d+l 



(d+l) d+ V Tmt [n 7 (T A )T A ]' 



(3.14) 



This choice of J(N) means that Ei nt cuts off rapidly for 
small chains, as it should, and that the contribution of 
the chains to the activation energy vanishes smoothly at 
T A - A simple way in which to interpolate between the 
small- N and large- N limits is to write 



J(N) 



:N d+1 



(3.15) 



This interpolation formula will be used extensively in the 
following analyses of experimental data. 

It seems natural to assume that £ R* (T) sets the length 
scale for the heterogeneous pattern of glassy domains and 
mobile regions discussed in SecHU The XC picture im- 
plies that molecular configurations confined within do- 
mains smaller than £ R* (T) are frozen because excitation 
chains are always subcritical at those length scales. On 
the other hand, domains appreciably larger than £ R*(T) 
can support XC-induced rearrangements and thus may 
be unstable against breaking up into smaller domains. I 
therefore propose that the R* (T) emerging from the XC 
theory be the same as the the R* (T) in the CTRW anal- 
ysis up to a dimensionless proportionality factor that, for 
present purposes, I set to unity. That is, I assume that 
this proportionality factor is absorbed into other quan- 
tities such as £ and the jump-length parameter a intro- 
duced in Eq.([23|). 

In [2], |3j , I argued that the diverging length scale R* (T) 
implies that the configurational entropy extrapolates to- 
ward zero at T = To, and does this in a way that is 
consistent with the apparent relation between the jump 
in the specific heat and the fragility near this point. This 
thermodynamic part of the theory is interesting; but it 
is even more speculative than the dynamic part, and is 
not directly relevant to the main issues being addressed 
here. Nevertheless, it does provide a clue about how 
to evaluate the glassy and mobile fractions, Vg{T) and 
V M (T) = 1 - V G (T). My guess in @, i| was that the 
excess configurational entropy might be proportional to 
the mobile fraction of molecules Vm(T), which in turn 
might be proportional to the surface-to volume ratio for 
glassy domains. More specifically, I proposed that, for 
d = 3, the excess entropy be proportional to 3 h* B /R*(T), 
where the length Kb — £h* B is the thickness of the inter- 
domain boundary. To extrapolate this geometric relation 
to T = Ta, where R* and Vg vanish, I propose that 



V G {T) 



R*(T) 



R*(T) 



(3.16) 



As will be seen, this simple extrapolation seems to work 
fairly well. 
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FIG. 1: (Color online) Intermediate scattering functions for 
T — 293 if, k — 2 A' 1 (red circles) and 1.2 A -1 (blue tri- 
angles). The data points are taken directly from the graphs 
shown in Refs. [22l |23|]. The solid curves are the theoretical 
results described in the text. 



IV. LESSONS FROM ORTHO-TERPHENYL 

The most extensively studied laborator y gl ass is ortho- 
terphenyl (OTP), for which the viscosity [19l.l2Cl] and the 
diffusion constant plj are known down to about the glass 
temperature T g = 240 K. The self intermediate scat- 
tering function F s (k,t) for this material is known from 
neutron scattering measurements down to T — 293 K 
[22I . I23L I24] ] , i.e. down to about its mode-coupling tem- 
perature Tmc HEHfl- In the following paragraphs, I use 
these three sources of information as guides for estimat- 
ing the model parameters introduced in the preceding 
Sections. Although I do not refer explicitly to them, my 
analysis is informed by numerical simulations, especially 
those reported in [27j ■ 



A. Scattering Data 

Start by looking at the scattering data, which pro- 
vides preliminary information about the time and length 
scales. A selection from the data of [22|, [23[ is shown 
in FigQ] along with theoretical curves computed using 
Eqs. (|2.19l - 12.22[) . As seen in the figure, the decoupling 
between the fast (j3) and slow (a) modes seems clear at 
the largest scattering wave vector, k — 2A~ 1 , but is not 
sharp at k = 1.2 A -1 or, as will be seen, at higher tem- 
peratures. By necessity, I assume that T = T s = 293 K 
is far enough into the thermally activated, heterogeneous 
regime for the preceding CTRW analysis to be valid, at 
least at long enough times. 

Eas. (j2~T5l) and (j2"Tg)) imply that, for large A, 
nM{k* , u) effectively has a pole on the negative real axis 



7 




FIG. 2: (Color online) Diffusion Constant L>oo(T). The data 
points (red circles) are taken directly from the graphs shown 
in Ref. [2lj]. The solid curve is the theoretical fit obtained 
from the CTRW and XC theories. The dashed curve is the 
mode-coupling approximation. 



at 



-A [1- P W /(A*)] 



-w . 



(4.1) 



At that point, ipoiu) is negligibly small but has a nonzero 
imaginary part, so that the integrand obtained by closing 
the contour around the cut on the negative real u axis has 
an approximate (^-function at u — —wq. This piece of the 
inverse Laplace transform contributes a rapidly decaying 
- i.e. "/?" - part of h M (k*,t*): 



n 



M ) 



A M e 



-w t* 



A 



M 



1 - V M f{k* 



(4.2) 



The function fi(j(k*,u) also has a pole at — wq; but its 
residue is vanishingly small for large A because molecules 
frozen in glassy domains do not undergo prompt /3 relax- 



ation. The plateau in F s (k*,t*) starts at t* 



+* 



1/wq 



and extends to t* = t* ~ 1. The height of the plateau is 



F? at =V G + {l-A M )V M = 



V G 



1 - V M f(k* 



(4.3) 



The OTP scattering data for T = T S = 293 A shown in 
Fig. Q] suggest that Ff at = 0.4 for k = 2 A~ x and (with 
much less certainty) Ff lat = 0.6 for k = 1.2 A -1 . Insert- 
ing these estimates into Eo. ([4.3p , I find Vm{T s ) — 0.75 
and aR*(T s ) = 0.60. As explained earlier, we cannot 
expect this CTRW analysis to describe /3 relaxation ac- 
curately at small times, and, indeed, the data in Fig. [1] 
do not exhibit a sharp fall-off at tp. The best I can do 
is to estimate tp(T s ) ~ 10~ 12 s for k = 2 A^ 1 . Knowing 
Vm and f(k*) at this point, I deduce from Eq. (|4.ip that 



FIG. 3: (Color online) Newtonian viscosity r?jy (T) . The data 
points (red circles) are those given in Refs. |T^, [20] . The solid 
curve is the theoretical fit obtained from the STZ and XC the- 
ories. The dashed curve is the mode-coupling approximation. 



tm ~ 2 x 10 14 s. Then, if tm is a thermal activation 
time of the form 



tm = to exp(T M /T), 



(4.4) 



and if t = 10~ 15 s, it follows that T M ~ 500 K. Thus, 
T = 293 K is inside the region of thermally activated 
dynamics, but only marginally so. 

The scattering data at T s = 293 K, k = 2 A -1 , also 
imply very roughly that t* a /t% = wqt^ ~ 10 3 . Again, 

knowing Vm and f(k*), I deduce that t a /tp = T a /rM — 
exp[T z /T s + a(T s )]/exp(T M /T s ) - 10 3 , where a(T s ) is 
defined in Eq.([Sl]). Thus, if T M = 500 if, then T z = 
2500 AT — T s a(T s ), which is substantially smaller than 
the value T z = 3100 A obtained in [H| by fitting the 
the viscosity data directly to an expression of the form 
of Eq. (|3.9p - a difference that turns out to be significant. 



B. Diffusion and Viscosity Coefficients 

The next step is to find model parameters that fit the 
diffusion and viscosity data shown in Figs. and [3] re- 
spectively. To evaluate the diffusion constant, start with 
the general formula for the mean-square displacement of 
a tagged molecule is 



'(**)> 



d 2 F s (k\t* 
dk* 2 



(4.5) 



where F s (k*,t*) is defined in Eq. (|2.20[) . To evaluate this 
formula in the limit of large t* , integrate around a small 
circle at the origin in the u plane, and use the expansion 



8 



i/jg( u ) ~ 1 — 3it/2 + • • •. The resulting integrand has a 
second-order pole at the origin, implying that 



(r* 2 (0> 



3aH 



(3/2)7> G + l/A" 



(4.6) 



Returning to dimensional units using Eqs. (|2.ip and (12. 4p . 
and defining the asymptotic diffusion constant by 



Doo i, we find 



£>oo(T) = 



3 a 2 £ 2 



(3/2)V G (T) Ta {T) + TM (T) 



(4.7) 



Thus, Z?oo oc £ 2 /r a at low temperatures where Vq — 
1 and r a ^> tm- When T approaches or exceeds T^, 
however, Vg vanishes and oc £ 2 /tm- 

An analogous formula for the Newtonian viscosity 
Vn(T) has been derived from the shear-transformation- 
zone (STZ) theory in 29]. It has the form 

m{T)=m^- [r' a (T)+T E (T)]. (4.8) 

This equation is the same as Eq.(7.10) in [29[ except that 
the latter contains an incorrect and almost entirely ir- 
relevant factor of 1/2. In Eq. (14.8p . t?o is a dimensional 
prefactor, r' a (T) is the same as r Q (T) in Eq. (|3.9p with T z 
replaced by T' z \ and 



te{T) = r exp 



rpl 
1 Z 



T 



(4.9) 



The quantities ks T' z and ks Te are, respectively, the 
characteristic formation energy for STZ's and the Eyring 
activation energy for STZ transitions. 



C. Stokes-Einstein Violation 

The results of Mapes et al.[2l[ indicate a strong vio- 
lation of the Stokes-Einstein relation; that is, the mea- 
sured value of the product D X (T) r)jy(T)/T is larger by 
a factor of about 10 2 near T g than it is near and above 
Ta (or Ta/c)j where it does remain constant as a func- 
tion of T as it should in a liquidlike material. The con- 
ventional explanation for this behavior is that the low- 
temperature heterogeneities allow a diffusing molecule to 
move further than it would in a homogeneous system. 
Such an enhancement of the diffusion length is built into 
the CTRW analysis in SecCH e.g. in Eq.(23]). How- 
ever, the effect is exactly cancelled out in Eq. (|4.7p by 
the correspondingly longer times that a molecule spends 
in glassy domains between diffusive jumps. There is, of 
course, no reason to believe that the Stokes-Einstein re- 
lation applies to the super-Arrhenius behavior of glass 
forming materials, where the mechanisms producing vis- 
cous deformation (localized shear transformations) and 
diffusion (dclocalized hopping) are quite different from 
one another. 



The preceding discussion of the characteristic temper- 
ature Tz suggests that the answer to the Stokes-Einstein 
question is that the Tz appearing in the diffusion analysis 
leading to Eq. (|4.7p is substantially smaller than the T' z 
appearing in the viscosity formula, Eq. (|4.8|) . In other 
words, the formation energies for density fluctuations 
may be substantially smaller than those required for cre- 
ating STZ's, which apparently are larger objects. (See 
the discussion in [2{|, Sec. IV.) The excitation-chain anal- 
ysis in Sec lIIII implies that the super-Arrhenius function 
a(T) and its ingredients are the same for both diffusion 
and viscosity; but there is no implication that the forma- 
tion energies Tz and T' z should be identical. Note then 
that the estimate 



Tg " T a (T g ) 



exp 



T' 
1 z 



To 



T 



10 2 



(4-10) 

implies that T' z — Tz~ 1000 if, in rough agreement with 
the estimate made by comparing time scales for decou- 
pling at T w T M c- 



D. Parameters and Length scales 

These order-of-magnitude estimates provide a starting 
point for using the diffusion and viscosity data to de- 
termine the remaining model parameters. Well within 
the super-Arrhenius region, the formulas for (T) and 
tjn(T), in Eqs. fl4.7j) and (|4.8p respectively, are most sen- 
sitive to the function a(T) which, in turn, is determined 
primarily by the choice of the Kauzmann/Vogel-Fulchcr 
temperature To and the product Ti nt 7q ■ This sensitiv- 
ity is best seen in the approximate formula for a(T), 
Eq. (|3.1 2p : To moves the Vogel-Fulcher singularity and 
Tint 7o determines its strength. I find that suitable val- 
ues of Ti n t 7o are always in the neighborhood of 10 s if. 
In order that Tint have roughly the same size as the for- 
mation energies T z and T' z , we need 70 = 3. I also find 
that T 178 K. 

Accordingly, the theoretical curves in Figs. and [3] 
have been plotted using T = 178 K, T in tJo = 10 5 , and 
70 = 3. Other parameters are: d — 3, v — ln(6), Ta = 
340 if, T M = 500 K, T E = 100 K, tj = 6 x I0 9 Pa s, 
T z = 2000.fr (for the diffusion constant), T' z = 3000 if 
(for the viscosity). With these parameters, A(T) > 100 
for all T < Ta, consistent with the large- A assumption 
made in earlier estimates. From the discussion following 
Eq. (|4.3p . I find that the parameter determining the jump- 
length distribution in Eq. (|2~5]l is a = 0.074. The fit to 
the diffusion constant requires a£ — 0.05 A; therefore I = 
0.7 A. Finally, knowing that V G (T S = 293 if) = 0.25, 
I find from Eq. (|3~TB]) that h B = 4.7 1 3.3 A. Fig. g| 
shows a graph of £ R* (T) in Angstroms as determined by 
these parameters; the corresponding graph of Vq(T) is 
shown in Fig. [5] 

Note that hg is approximately half the diameter of 
an ortho-terphenyl molecule. However, £ is only about 
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FIG. 4: Length scale £R*(T), in Angstroms, as a function of 
temperature T, for the OTP parameters determined in the 
text. 



a tenth of this diameter, implying that the elementary 
links in excitation chains and the individual diffusive 
jumps in mobile regions are quite small, and therefore 
that those motions might involve cooperative rearrange- 
ments of many molecules rather than being simple, pair- 
wise exchanges of positions. Interestingly, this value of 
I is of the order of the Lindemann length, and perhaps 
is consistent with the idea that the excitation chains de- 
scribe a local melting of the system during configura- 
tional rearrangements. Other authors such as H3| 
have deduced similarly small length scales. For example, 
using the parameters listed above and T g = 240 K, I find 
that £R*(T g ) = l&Al = 111, or less than 2 molecular 
diameters, roughly consistent with the results of [3(3] (de- 
pending on - among many other uncertainties - whether 
R* is interpreted as the diameter or radius of a glassy 
domain) . 

The theoretical fits to experimental data shown in 
these figures are interestingly imperfect. In some cases, 
they can be improved cosmetically by adjusting parame- 
ters such as Tm, Tz, etc., or by arbitrarily adjusting val- 
ues of Vg for T > 293 K slightly downward from those 
determined by Eq. p.l6[) . Moreover, if I relax the con- 
straint that the product Ti„ t 7q be the same for both 
diffusion and viscosity, and instead let this quantity be 
0.8 x 10 5 for diffusion and 1.2 x 10 5 for viscosity, the 
theoretical curves for both D^T) and 1]n(T) become 
indistinguishable from the data over their entire tem- 
perature ranges. I doubt that such an adjustment is 
physically realistic. To emphasize this point, in Figs. 
[5] and [31 I have added mode-coupling approximations of 
the form oc (T - T M c) 2 and r) N cx (T - T M c)~ 2 , 
with Tmc = 290 K. These simple approximations fit the 
high-temperature data extremely well, and illustrate the 
way in which the experimental results cross over at about 



FIG. 5: Glassy probability Vg{T), as a function of tempera- 
ture T, for the OTP parameters determined in the text. 



Tmc from solidlike, thermally activated behavior at low 
T to liquidlike behavior at high T. 

The discrepancies in fitting the scattering data are in 
part due to the same difficulty - that the crossover re- 
gion at temperatures near Tmc seems theoretically in- 
tractible. To illustrate this situation, I show in Fig. [6] 
the scattering functions F s (k,t) for k = 2 A~ x and tem- 
peratures T = 327 K, 306 K, 293 K (with accompanying 
experimental data) and T = 280 K (for which no data 
are available). The two higher temperatures are clearly 
beyond the range of validity of the activation theory; the 
CTRW analysis is too blunt an instrument to resolve the 
details of the (3 relaxation, and the activation barriers are 
too low. Nevertheless, the agreement with experiment, 
even at T = 327 K, seems qualitatively reasonable, im- 
plying that the temperature and time scales may be at 
least approximately correct. 



V. STRETCHED EXPONENTIAL 
RELAXATION 

Finally, it is useful to look more closely at the long-time 
tails of the scattering functions shown in Figs. Q]andE] 
My interest here is not so much in the experimental data 
for OTP as it is in the implications of the theory for in- 
terpreting extensions of that data. I find that the CTRW 
theory predicts stretched-exponential decay of F s (k*,t*) 
with a temperature and wavenumber dependent index b 
that crosses over to the intrinsic b = 0.5 only at unob- 
servably long times t*. 

Especially in the case of the curve for T = 293 K, 
k = 2 A~ X , the CTRW analysis should be accurate at 
large t*. I have replotted that curve in Fig. [71 here in 
the form — log 10 [— log 10 F s (k, t)] as a function of log 10 (f), 
so that the slope is asymptotically equal to (minus) the 
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"og 10 (t) (s) log 10 (t) (s) 



FIG. 6: (Color online) Intermediate scattering functions for 
k = 2 A' 1 , T = 327 K, 306 if, 293 K, and 280 if (dashed 
line), from left to right. 

strctchcd-exponential index b. For comparison, I also 
have plotted the analogous curve for T = 327 K, again 
with k — 2A~ 1 , although that temperature is well into 
the mode-coupling region according to Figs. [2] and [3] 
The experimental points are shown for both cases. 

These curves, and others like them (not shown here) for 
different values of k and T, exhibit many of the features 
seen in Fig. 2 of Ref [l|. (The dimensionless wavenumber 
k* was denoted simply by k in that paper.) The analysis 
in [l[ tells us that, in the limit of small k* at fixed large 
t*, the behavior becomes diffusive, b — > 1. Conversely, at 
fixed nonzero k*, in the limit of large £*,£>—> 0.5. Fig. 
2 of [l[ also shows that the apparent value of b changes 
continuously from b = 0.5 at large k* to b — 1 at small 
k*. 

These observations provide a clue about what to ex- 
pect in the temperature dependent theory. Because 
k* = kt R* (T) , the wavenumber k* in [lj is a proxy for 
the temperature. At fixed scattering wavenumber k, k* 
decreases as T increases. If b approaches unity as k* de- 
creases toward zero in the T-independent theory, then we 
may expect that, at fixed k, b will behave in the same way 
with increasing T as R*(T) becomes small. Apparently 
this is what happens; but the temperature dependence 
of b turns out to be very slow at low T, and the inter- 
esting high-temperature behavior is beyond the limits of 
validity of the present theory. 

I have pursued this analysis in Fig. [7]by comparing the 
CTRW results with pure stretched-exponential functions 
of the form 

F s SE (k*,t*) =FP lat exp(-C SE t* b ), (5.1) 

where Cse and b are fitting parameters and Ff lat is de- 
termined by Eq. (|4.3j) . For the upper curve, T — 293 K, 
Cse = 0.75, and b = 0.70. For the lower curve, 



FIG. 7: (Color online) Double logarithmic plot of the interme- 
diate scattering function for k — 2 A" 1 , for T — 293 K (upper 
black curve) and T = 327 K (lower black curve) . The red and 
blue dashed curves are fits to pure stretched exponentials for 
the upper and lower curves with indices b = 0.69 and 0.73 
respectively. The red circles and blue triangles are the same 
data points as those shown in Fig. [6] 



T = 327 K, C SE = 0.65 and b = 0.72. The fits seem 
to be exact all the way from the beginning of the plateau 
out to the point where the slope changes to b = 0.5 sev- 
eral decades in time later. At this crossover, in both 
cases, F s has dropped from about 0.4 to an unobserv- 
ably small value of the order of 10 -5 . To see the intrinsic 
b = 0.5, or the value of b appearing in a generalized 
form of the waiting-time distribution ipo(t*) discussed 
following Eq. (|2.2[) , we would have to go well beyond this 
crossover point. Thus, any attempt to deduce b from 
an experimental curve like this would produce an inter- 
mediate value, equal in this case neither to unity nor to 
0.5. 

Further comparisons indicate that the effective index 
b deduced in this way increases slowly with increasing T. 
The value of b at T = 327 K is slightly larger than that 
at T = 293 K. At the higher T, b is measured over a 
smaller range of times and therefore with less precision; 
but the small difference between the b values is within the 
accuracy of the fits. However, at a lower temperature, 
T = 280 AT, the parameters b and Cse are the same as 
those at 293 K to the accuracy of my calculation, b does 
depend more strongly on k at fixed T. At k = 1.2 A -1 
and T = 293 K, I find b = 0.77 - appreciably larger 
than at k — 2A~ 1 . I also have checked to make sure 
that the value of b measured by this fitting procedure is 
always unambiguously equal to unity if i/jcit*) is a simple 
exponential function. 

To see what is happening here in somewhat more de- 
tail, look at an approximation for F s (k*,t*) that is accu- 
rate for A(T) 3> 1 and for values of t* large enough that 
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we can use the small-iu approximation, ^q(—w ± ie) « 
l+3w/2+- • B(uj), where B{w) is defined in Eq. ([2T2T]) . 
I find 

w A (k*,T) f°° dw 



F a (k*t*) = 

/(**) Jo irw 
(2/3) exp(-wi') 

X [ w -^(fc«,r)] 2 + [(2/pH] ; 



where 



w A (k*,T) 



2 [1 - f(k*)\ 
3V G (T)f(k* 



(5.2) 



(5.3) 



Apart from the prefactor /(fc*) -1 , which does not affect 
6, all of the temperature and wavenumber dependence 
is contained in WA(k* ,T). That is, the effective b(k,T) 
is a function of only wa within the range of validity of 
Eq.UEU)- 

In the limit of vanishingly small k* , and at large fixed 
t*, f(k*) « 1, B{w A ) < 1, and the integrand in Eq.(|B"^]) 
is accurately approximated by a 5 function at k; = oc 



Therefore 



F s (k*,t*) 



-W A t 



exp 



2a 2 e k 2 t 



3V G (T) Ta (T) 



, (5.4) 



which confirms that b — > 1 in the small-fc* limit. 

At fc = 2 A" 1 and T = 293 K : however, w A = 2.8; 
therefore, for times t* of the order of unity or greater, 
the integrand in Eq. l|5.2p is negligibly small at w ~ wa- 
At values of t* so large that the integrand is dominated 
by the product B(w) exp(—wt*) oc exp(— wt* — 1/w), 
there is a sharp peak at w = 1/ Vi* , and a saddle-point 
approximation yields 



F s {k*,t*) oc2Vt 



,-2 \/F 



(5.5) 



At large £*, this function is essentially the same as 4>a(t*) 
given in Eq. (|2.8|) . However, the approximations made in 
obtaining this result are not valid until the saddle point is 
well below the maximum of B(w) at w — 2/3. Therefore, 
the saddle must occur at values of w less than about 0.3, 
and Ea.(|5.5p is correct only for t* > 10, i.e. t > 10~ 7 
s, which is roughly where the crossover from b = 0.70 to 
b = 0.5 occurs on the upper curve in Fig. 

The remarkable feature of these numerical results is 
that the scattering function F s (k*,t*) predicted by the 
CTRW theory is so well fit by a simple stretched- 
exponential in its experimentally observable domain. I 
have checked that Eq. (|5.2[) accurately reproduces the 
large-i* behavior of F s (k*,t*) in both the effective b(k, T) 
regime and the asymptotic b — 0.5 limit, but can find 
no hint from this formula that Eq. (|5.ip is anything 
more than a simple - and therefore non-unique - two- 
parameter fit to a smooth function that decreases more 
slowly than exponentially in a restricted range of values 
of t*. The analysis does make it clear that stretched- 
exponential behavior depends on there being an intrinsic 



stretched-exponential waiting time distribution tpc(t*). 
The important outstanding question is whether the in- 
trinsic b itself must be temperature dependent in order 
that the observed index b exhibit physically realistic tem- 
perature and wavenumber dependence. 



VI. CONCLUDING REMARKS 

The introduction to this paper begins with a list of 
questions about the role of spatial heterogeneities in the 
dynamics of glass- forming liquids. Several of those ques- 
tions were addressed and at least partially answered in 
the preceding paper [lj; others are discussed in this one. I 
conclude here by summarizing my present opinions about 
each of them. 

Structural or dynamic heterogeneity? I think that this 
is largely a semantic issue. We referred in [l[ to papers 
by Widmer-Cooper and Harrowell 3l| , who showed that 
molecules in mobile regions - regions of high "propensity" 
- have smaller than average Debye- Waller factors, mean- 
ing that they are participants in elastically soft modes. 
Such modes are ultimately structural in nature, but are 
to some degree delocalized and not easily detected by 
looking at near-neighbor structural correlations. I think 
that the mobile regions are the locations of these soft 
modes, and that the nonlinear /3 rearrangements occur 
predominantly in these places where the linear elastic 
modes are softest and most highly excited by thermal 
fluctuations. The fact that the pattern of mobile regions 
and glassy domains is fluctuating on a time scales adds 
an extra dynamic element to the picture of the heteroge- 
neous system as a whole. 

Heterogeneity and non-Gaussian displacement distri- 
butions? The conclusion reached in |l|, consistent with 
the results of |l3j, is that non-Gaussian displacement 
distributions are robust, mathematical features of essen- 
tially all heterogeneous diffusion models in which there is 
strong decoupling between the fast and slow modes. This 
result depends only on the requirement that ipG(t*) de- 
cays rapidly enough that ipG(u) is very small at large u; 
it does not depend at all on whether ipait*) is a stretched 
exponential function. I have not repeated this analysis 
here because, apart from some inconsequential technical 
differences, it is exactly the same as in Sec. V of [l[. Note 
that, while stretched exponentials and non-Gaussian dis- 
tributions are both generated by spatial heterogeneities, 
these are independent phenomena; non-Gaussian diffu- 
sion can occur in the absence of stretched exponential 
relaxation. 

Violation of Stokes- Einstein? Although it might have 
been satisfying to find a theoretically more elegant expla- 
nation - longer diffusion paths or a breakdown of ergodic- 
ity - 1 argue in Sec. IIVI -C that the violation of the Stokes- 
Einstein relation has a more prosaic origin. Specifically, 
as a glass-forming liquid becomes more solidlike below 
Ta, the localized shear transformations producing vis- 
cous deformation become less and less like the delocalized 



12 



hopping mechanisms responsible for diffusion. Looked at 
from this point of view, the large violation of Stokes- 
Einstein is a prediction of the shear-transformation-zone 
(STZ) theory [IH, and the scale of the low-T violation 
(the factor of 100 found in [2l| ) is independently pre- 
dicted by the high-T time-scale analysis described in Sec. 
IMA. 

Length scales? I cannot prove that the estimates of 
relatively small length scales in Sec. IIVI D are unique 
interpretations of the available data, but I suspect that 
the conclusions are at least qualitatively correct. If so, 
it seems especially urgent to understand what is happen- 
ing here. What kinds of molecular rearrangements cor- 
respond to such small net displacements? On the other 
hand, there are still substantial uncertainties. For exam- 
ple, the poorly understood parameter 70 plays an impor- 
tant role in these estimates, and the quantity R*(T) that 
appears as the critical size of an excitation chain in Sec. 
IIIII may not be exactly the same as the size of the glassy 
domains in Sec. [TXJ The most serious experimental un- 
certainties pertain to the scattering data. It would be 
useful if these could be extended to lower temperatures 
and longer times. 

Heterogeneity and stretched exponentials? The connec- 
tion between heterogeneity and stretched exponentials 



presented in [l[ seems plausible; but it is only one step 
more specific than the conventional picture in which a 
stretched-exponential function is just an incoherent su- 
perposition of ordinary exponential decays occurring in 
different parts of a heterogeneous system. I think it may 
be possible to make a first-principles estimate of the in- 
trinsic, possibly T-dependent, stretched-exponential in- 
dex b that appears in the glassy, waiting-time distribu- 
tion ipG{t*)- However, the discussion in Sec. fVl implies 
that this intrinsic index is impossible to measure directly, 
and that further calculations are needed in order to make 
contact with experiments. The positive aspect of the sit- 
uation is that the theoretical correlation functions, like 
many seen in experiments and numerical simulations, ex- 
hibit sharply defined, temperature and wavenumber de- 
pendent indices b(k, T) over long times. Those indices 
should be calculable from the intrinsic index b via CTRW 
analyses such as the one described here. 
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